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ABSTRACT 

We describe new radial velocity and X-ray observations of extremely low- 
mass white dwarfs (ELM WDs, ~ O.2M ) in the Sloan Digital Sky Survey Data 
Release 4 and the MMT Hypervelocity Star survey. We identify four new short 
period binaries, including two merger systems. These observations bring the to- 
tal number of short period binary systems identified in our survey to 20. No 
main-sequence or neutron star companions are visible in the available optical 
photometry, radio, and X-ray data. Thus, the companions are most likely WDs. 
Twelve of these systems will merge within a Hubble time due to gravitational 
wave radiation. We have now tripled the number of known merging WD sys- 
tems. We discuss the characteristics of this merger sample and potential links to 
underluminous supernovae, extreme helium stars, AM CVn systems, and other 
merger products. We provide new observational tests of the WD mass-period 
distribution and cooling models for ELM WDs. We also find evidence for a new 
formation channel for single low-mass WDs through binary mergers of two lower 
mass objects. 
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INTRODUCTION 



Before the re cent discovery of th e six short period binary WD systems containing 



0.2 M ? ELM WDs JKilic et al.lhoiOblfllMullallv et al.l2009l:lMarsh et al.l boiol;lKulka rni fc van Kerkwiik 



2010; Badcnes et al 



2009 



Steinfadt et al. 



2010; 



Kawka et al. there were only six dou- 



ble WD systems known to have short enough orbital periods to merge within a Hubble time 



sec 



Nelemans et al.ll2005l . and references therein). Even though the Supernovae la Progen- 
itor Survey (SPY) identified f»100 new double degenerate systems from a sample of 1000 



WDs observed using the Very Large Telescope, iNapiwotzki et al.l ( 120041 ) identify only eight 
merger candidates (including both WD+WD and sdB+WD systems) from the SPY and 
the literature. This fraction (8 out of 1000) is relatively small co mpared to the fraction of 



merge r systems (6 out of 8) among the ELM WDs in the literature (jKilic et al.ll2007bl . 12009 
2010bl Jcl). 



Low-mass (M < 0.45M©) He-core WDs are formed when a companion strips the outer 
envelope from a post main-sequence star before the star reaches the tip of the red giant 
branch and ignite s helium. T hey are usually found in close binaries, mostly double de- 
generate systems ( lMarshlll995l ). However, about half of the low- mass WDs in the field do 



not show any radial ve lo city variations, ind icating that they are single (IMaxted et al.l 12000 



Napiwotzki et al. 


2007 


)• 


Kilic et al. 


(2007c) 



from old metal-rich stars that truncate their evolution prior to the helium flash from severe 
mass loss. They estimate a binary fraction of 50% for ~ 0.4 M Q WDs. However, they predict 
that the binary fraction rises to 100% for ~ O.2M ELM WDs, since such extreme mass loss 
rates are not expected even for the most metal-rich stars in the Galaxy. 



1.1. The ELM Survey 



To study the binary fraction of ELM WDs and to find future merger systems, we have 
undertaken an MMT survey of all previou sly identified ELM WDs from the SP SS DR4 area 



( Eisenstein et al.ll2006t iLiebert et al.ll2004l ) and the Hypervelocity Star survey oflBrown et al. 



( 120061 ). The latter sample is described in a com panion paper by IBrown et al.l (j2010al ). Part 
of the SDSS DR4 sample was discussed earlier in lKilic et al.l (j2010bf ). Here, we present radial 
velocity observations of 7 more ELM WD candidates from the SDSS sample. The combined 
sample contains 23 WDs. We also report on Chandra X-ray observations of three ELM WDs 
to search for X-ray evidence of milli-second pulsar (MSP) companions. 



While almost all SDSS ELM WDs show radial velocity variatio ns, none show evidenc e 
of main-sequence or neutron star companions in the optical, radio ( lAgiieros et al.l l2009al ) . 
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or X-ray (this paper and lAgiieros et al.l l2009bl ). Hence, they are almost certainly binary 
WD systems. Out of the 20 short period binary systems in our sample, 12 have merger 
times shorter than a Hubble time. Our ELM survey has now tripled the number of known 
merging WD systems. Based on our sample and previously known short period binary WD 
systems, we study the mass-period distribution and future evolutionary prospects, including 
the connections to unusual stellar populations and underluminous supernovae (SNe). 

Our paper is organized as follows. In Section 2 we discuss our radial velocity and X-ray 
observations. In Section 3 we present an analysis of the optical spectroscopy and X-ray data 
and the nature of the companions. In Section 4 we discuss the sample characteristics of 
ELM WDs, including the binary fraction, mass-period distribution, merging systems, and 
the implications for underluminous SNe and other unusual merger products. We conclude 
in Section 5. 



OBSERVATIONS 



2.1. Optical Spectroscopy 



Eisenstein et al.l (120061 ) identified 13 low-mass WD candidates with log g < 7. iKilic et al. 
(j2010bl ) discuss follow-up observations of four of these targets. We selected 6 more targets for 
follow-up MMT spectroscopy. These targets are SDSS J002228.45+003115.5, J002207.65-101423.5, 
J123410.36-022802.8, J162542. 10+363219.1, J204949.78+000547.3, and J225242.25-005626.6. 
We also obse r ved S PSS J234536.48-010204.8, another ELM WD candidate identified by 
Liebert et all (booj ). 



There are three more targets in the lEisenstein et al.l (120061 ) ELM WD sample that lack 
follow-up spectroscopy. These targets (SP SS J1056+6536, J1426+ 0100, and J1630+4233) 
have log g > 6.9 based on the model fits by lEisenstein et al.l (120061 ) and log g ~ 7 based on 
the fits by IKilic et al.l (j2007al ). Radial velocity observations of these objects will be useful 
for the identification of more short period binary WP systems. 

We used the 6.5m MMT equipped with the Blue Channel Spectrograph to obtain mod- 
erate resolution spectroscopy of seven ELM WP candidates over 17 different nights between 
2008 September and 2010 July. We operate the spectrograph with the 832 line mm -1 grating 
in second order, providing wavelength coverage 3600 A to 4500 A and a spectral resolution 
of 1.0 - 1.2 A, depending on the slit size used. We obtain all observations at the parallactic 
angle, with a comparison lamp exposure paired with every observation. We flux-calibrate 
using blue spectrophotometric standards (IMassev et al.lll988f). and we measure radial veloc- 
ities using the cross-correlation package RVSAO (IKurtz & Minklll998l ). Our radial velocity 
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measurement procedures are described in iKilic et al.l (120091 l2010bf ) 



2.2. X-ray observations 

2.2.1. Motivation 

X-ray emission from radio MSPs is expected to be less highly beamed than radio emis- 



Zavlin et al. 


[2002; 


Zavlin 


2006) 



companion to the ELM WDs will be gravitationally bent, allowing observation of >75% of 
the neutron star surface in X-r ays even if the rad io pulsar beam (produced well above the 



surface) misses our line of sight (lBeloborodovll2002l ) . All 15 radio MSPs with precise positions 



i n unconfused regions of the globular clu ster 47 Tuc have been clearly detected in X-rays 



(IHeinke et al.ll2005l ; iBogdanov et al.ll2006l ). with no observed correlation between their radio 
and X-ray luminosities (expected, as they are produced via different mechanisms). There is 
also no evidence for dif ferences in the X-ray properties of globular cluster vs. Galactic MSPs 



(IBogdanov et al.ll2006f ). 



This result allows us to use the 47 Tuc MSP sample (with accurate Lx values due to 
its well-known distance) to predict that other MSPs should have X-ray luminosities above 
Lx (0.5-6 keV)= 2 x 10 30 erg/s, the minimum Lx of MSPs in 47 Tuc. Thus, deep X-ray 
upper limits can rule out the presence of an MSP, as in the ELM WD SDSS J0917+4638 
(jAgiieros et al.l l2009bl ). The MSPs in 47 Tuc show X-ray spectra dominated by thermal 
bla ckbody-like emission fr om the neutron star polar cap surfaces, with temperature 1—3 x 10 6 
K (IBogdanov et al.l 120061 ) . Some MSPs also show additional (sometimes much higher Lx) 
components due to magnetospheric emission or pulsar wind shocks, but there is no evidence 
for a MSP lacking the thermal component. 

We obtained Chandra observations of three ELM WDs (SDSS J105353.89+520031.0, 
J123410.36-022802.8, and J092345.59+302805.0) to search for MSP X-ray emission. None 
of these WDs had been previously observed in the X-rays since the ROSAT All-Sky Survey 
( IVoges et al.l 119991 ) . where the short exposure time was insufficient to place useful limits. 
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2.2.2. Data Analysis 



We observed the three systems with Chandra 1 ?, ACIS-S instrument in Very Faint mode, 
for the exposure times noted in Table 1. We used CIAO 4.2 (with CALDB 4.2.1)0 to 
reprocess the observations with current calibrations and reduce the backgrounds using Very 
Faint mode cleaning. We constructed images in the 0.3-6 keV band, and found n o X-ray 
photo ns within the 1" error circles around each source. We use Poisson statistics (IGehrels 
19861 ) to place 99% confidence u pper limits o n the c ount rate. We compute distances to the 
ELM WDs using the models o f IPanei et al. ( 20071) and the SDSS g magnitudes, and Nh 
values using the Colden too@ (jDickev fc Lockman 1990 ). We use PIMM£§ and the X-ray 
spectrum of the faintest MSP in 47 Tuc (47 Tuc-T, 134 eV blackbody) to produce 0.5-6 keV 
Lx upper limits, which we list in Table 1. 

The 99% confidence Lx upper limits are all at least a factor of 10 fainter than the 
faintest MSP identified in 47 Tuc, and 21 times fainter than the median Lx of these MSPs 
(jBogdanov et al.l 120061 ). Thus the lack of X-ray emission from these WDs is strong evidence 
that the unseen companions are not MSPs. 



3. RESULTS 

Four of our seven targets exhibit significant radial velocity variations, with peak-to-peak 
velocity amplitudes between 124 and 232 km s _1 . We weight eac h velocity measurement b y 



its associated error and solve for the best-fit orbit using the code of iKenyon fc Garcia I (119861 ). 
Figures 1-6 show the observed radial velocities and the best fit orbits for our targets. The 
heliocentric radial velocities are best fit with circular orbits with periods 1.9-11.8 hr. We 
present the best-fit orbital period (P), semi-amplitude (K) of the radial velocity variations, 
systemic velocity (which includes a small gravitational redshift term), the time of spectro- 
scopic conjunction, and the mass function in Table 2. J2049+0005 is excluded from this list 
due to its low surface gravity, which implies that it is not a WD (see below). 

Two of the stars in our sample, J2252— 0056 and J2345— 0102, do not show significant 
velocity variations. The best-fit orbital solutions have semi-amplitudes 25-43 km s _1 . These 
two objects may be low inclination systems or just single stars. 



1 http://cxc. harvard . edu/ciao/ 
2 http://asc. harvard.edu/toolkit/colden.jsp 
3 http: / / asc.harvard.edu/toolkit / pimms.jsp 
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We perform model fits to each individual spectrum and also to the average composite 
spectra using synthetic DA WD spectra kindly provided by D. Koester. This model grid 
covers temperatures up to 30,000 K. We use an NLTE grid computed by I. Hubeny (see 
Allende Prieto et al.l 120091 ) to fit the spectrum of J2345— 0102, which requires a solution 
slightly above 30,000 K. We use the individual spectra to obtain a robust estimate of the 
errors in our analysis. Figure 7 shows the composite spectra and our fits using the entire 
wavelength range. The best-fit T eS and log g values are given in Table 3. We obtain best-fit 
solutions of 17, 890 — 33, 130 K and logg = 6.12 — 7.38 for six of our targets, confirming that 
they are low- mass WDs. 

Figure 8 shows the effective temperatures and surface gravi ties for our targets (red - 



circles) plus the previously identified ELM WDs in the lit erature (jKific et al.l l2007bl 12009 



Vennes et al. 


2009; 


Kilic et al. 


2010b: 


Brown et al. 


2010a) 



companions to milli-seco nd pulsars PSR J1012+5307 and J1911-5958A ( Ivan Kerkwijk et al. 



1996: 


Bassa et al. 


2006 


. see § 


WDs from 


Panei et al. 


(2007) 



2006 L see §4.4). Solid lines sh ow the constant m ass tracks for low mass 



our six WD targets have masses ranging from 0.20 to 0.42 M . 



3.1. Notes on Individual Objects 



3.1.1. J0022+0031 



The best-fit orbital period for J0022+0031 is 11.8 hr. However, the limited number of 
follow-up observations allow for an alias at 7.9 hr period. J0022+0031 is a binary system 
containi ng a 120 My r old .38 M & WD and a M > O.21M companion star. Based on the 
updated iPanei et al.l (120071 ) models, it has an absolute magnitude of M g = 9.8 mag and a 
distance of 790 pc. 

For all six WDs in our sample, we combine the spectra near maximum blue-shifted 
and red-shifted velocities into two composite spectra. If there is a contribution from a 
companion object, it may be visible as an asymmetry in the line profiles. We do not see any 
obvious asymmetries in the line profiles and our optical spectroscopy does not reveal any 
spectral features from companion objects. A M > 0.21 main -sequence star companion to 
J0022+0031 would have Mi < 10.0 mag fcroupa k Toutlll997t ). brighter than the low-mass 
WD [Mi ~ 10.4 mag) and detectable in the /-band. Such a main-sequence star companion 
is ruled out based on the SDSS photometry. 



Using the mean inclination angle for a random stellar sample, % = 60°, the companion 
has a mass of 0.26 M Q . There is a 2.5% probability that the companion is a 1.4-3 M & neutron 
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star. Given the small probability of a neutron star companion and the unsuccesful s earches 
for neutron stars around other ELM WDs (Section 2.2, see also lAgiieros et al.ll2009al lbl). the 
companion is most likely another WD. Short period binaries may lose angular momentum 
through gravitational wave radiation and merge within a Hubble time. The merger time for 
the J0022+0031 system is longer than a Hubble time. 



3.1.2. J 0022- 1014 

J0022— 1014 is a 1.9 hr period binary system with a 0.33 M Q WD primary and a M > 
O.19M secondary. The primary WD is 70 Myr old, it has M g = 9.2 mag, and is located at 
a distance of 1.3 kpc. A M > O.19M main-sequence star companion would have Mj < 10.3 
mag, 30% fainter than the ELM WD. Such a companion is ruled out based on the SDSS 
photometry. The companion to J0022— 1014 is most likely another WD, and specifically a 
low-mass WD. There is a 83% probability that the companion is less massive than 0.45 M Q , 
a He-core low-mass WD. This system will merge in less than 730 Myr. 



J0022— 1014 has a proper motion of fi a cos5 = —7.8 and us = —13.2 mas yr 1 fjMunn et al 



2004J). Based on the mass and radius estimates, the gravitational redshift of the WD is 
8.1 km s _1 . Therefore, the true systemic velocity is —46.6 km s " 1 . The ve l ocity com- 



ponents with respect to the local standard of rest as defined by iHogg et al.l ( 120051 ) are 



U = 95 ± 23, V = -53 ± 21, and W = 33 ± 8 km s _1 . J0022-1014 is a disk star. 



3.1.3. J 1234- 0228 



J1234— 0228 was originally identified as a very low mass WD by iLiebert et al.l (120041 ) 
and is included in the DR4 WD catalog of lEisenstein et al.l ( 120061 ) . Based on our MMT 
observations over 5 different nights separated by about a year, J1234— 0228 is a binary 
system with an orbital period of 2.2 hr. The visible component of the binary is a 70 Myr old 
WD with M g = 8.4 mag at a distance of 780 pc. For an inclination angle of 60°, the mass 
function requires a 0.11 M Q companion. A main-sequence companion of this mass would 
contribute 7% excess flux in the /—band. No excess flux is observed in the SDSS photometry 
and spectrum of th is object, therefore the companion is most likely another WD (see also 
Liebert et aDl2004h . 



There is a 94% chance that the companion is a low-mass (< 0.45M Q ) WD, and more 
interestingly there is a 85% probability that the combined mass of the two stars in this 
system is lower than 0.45 M Q . The merger time for the J1234— 0228 system is shorter than 
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2.7 Gyr. Therefore, within the next few Gyr, this system will merge and likely form a single 
low- mass WD. 



J1234— 0228 has a proper motion of /j, a cosS = —14.5 and \x& = —12.3 mas yr 1 ( IMunn et al 



20041 ). Based on the mass and radius estimates, the gravitational redshift of the WD is 
3.7 km s _1 . Therefore, the true systemic velocity is +46.6 km s " 1 . The ve l ocity com- 
ponents with respect to the local standard of rest as defined by [Hogg et al.l ( 120051 ) are 
U = -5 ± 14, V = -78 ± 12, and W = 21 ± 7 km s _1 . J1234-0228 is also a disk star. 



3.1.4. J 1625 +3632 

J1625+3632 has a best-fit orbital period of 5.6 hr. The MMT spectrum reveals a weak 
He I line at 4471 A. Our best -fit DA WD mode l has T cS = 23570 K and log g = 6.12, similar 



to the sdB star HD 188112 (IHeber et ajj 120031 ) . Based on these model fits, the primary is 



a ^160 Myr old 0.20 M Q WD at a distance of 2.8 kpc. J1625+3632 does not display a 
significant proper motion, but given its distance and systemic velocity, it is likely a member 
of the thick disk or halo. 

The relatively small semi-amplitude (K = 58.4 km s _1 ) of the radial velocity variations 
imply that the lower limit on the mass of the companion is 0.07 M & . A main-sequence 
companion of this mass would contribute less than 1% in the optical and would not be 
detected in our observations or the SDSS photometry. Given the accuracy of the SDSS 
photometry, any main-sequence companion more massive than about 0.17 M & would have 
been detected. Hence, only 0.07-0.17 M Q main-sequence companions are allowed based on 
the mass function and optical photometry. Near-infrared observations will be useful to rule 
out such companions. 

J 1625+3632 may be a low inclination system. However, the probability of a neutron 
star companion is only 0.6%. Given the narrow range of possible main-sequence companions 
and the unlikelihood of a neutron star companion, the companion is likely another WD. 
There is a 96% chance that the companion is less massive than 0.45 M , i.e., a low-mass 
He-core WD. 



3.1.5. J2049+0005 



Figure 9 shows the best-fit WD model to our MMT sp ectrum of J2049+0005. This star 
was identified as a T e g = 8660 K and lo g q = 5.48 WD by lEisenstein et al.l (120061 ). Further 
analysis of the same SDSS spectrum by lKilic et al.l (j2007al ) suggested that J2049+0005 has 
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log (7 < 5, and therefore is not a WD. Similarly, our model fits to the higher resolution and 
higher signal-to-noise ratio MMT spectrum find a best-fit solution that is exactly our lowest 
gravity model (log g = 5). This model poorly matches the high order Balmer lines, indicating 
that the surface gravity is lower than log g = 5 and that J2049+0005 is most likely an A 
star. 



3.1.6. J 2252- 0056 

Based on 30 different spectra, the best-fit orbital period is 10.3 hr with K — 25.1 ± 1.5 
km s _1 . However, the radial velocity variations are not significant. The best-fit orbit has 
a x 2 slightly better than a constant-velocity fit. If J2252— 0056 is a binary, it is likely 
a low inclination or a pole-on system. Given the 25 km s _1 upper limit to the velocity 
semi-amplitude, J2252— 0056 could be a binary with either month-long (assuming an 0.5 
Mq companion and i = 60°) orbital period or relatively pole-on inclination (M sin (i) = 
0.03 — O.O5M for 2-12 hr periods). There is no evidence of a companion in our MMT 
spectrum or the SDSS observations. 



3.1.7. J 23 4 5- 01 02 

J2345— 0102 is another object for which we detect no significant velocity variations. 
The best-fit orbital period is 7.2 hr with K = 43.0 ±0.7 km s , but the best-fit orbit is a 
slightly better fit to the data than a constant- velocity fit. Again, if J2345— 0102 is a binary, 
it may be a low inclination system. The observed 43 km s _1 upper limit to the velocity 
semi-amplitude corresponds to a binary with either month-long orbital period or relatively 
pole-on inclination (Msin (i) = 0.05 — 0.10M o for 2-12 hr periods). 

The MMT spectrum of this object is best-explained with a 23 My r old 0.42 M G) WD 



at a distance of 2.0 kpc. Based on 5 epochs from the USNO-B+SDSS, iMunn et al.l (120041 ) 
measure a proper motion of \i a cosb = —26.0 and fis = —36.6 mas yr _1 . Given its relatively 
large systemic velocity, proper motion, and distance, J2345— 0102 is clearly a halo star- 



About half of the known 0.4 M & WDs seem to be single (see the discussion in iKilic et al. 



2007d ). Hence, the lack of significant radial velocity variations is not very surprising for this 



system. 
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DISCUSSION 



4.1. The Binary Fraction of ELM WDs 



Kilic et al.l ( ]2007bl . l2010bl ). iBrown et al.l (j2010al ). and this paper discuss radial velocity 
observations of 11 ELM WD candidates identified in the SDSS DR4 (lEisenstein et al.ll2006l) 
and 12 ELM WDs identified in the Hypervelocity Star Survey of IBrown et al.l ( 120061 ). One of 
the targets is common t o both samples. I n addit i on, one of the DR4 target s (J2049+0005) 



i s mos t likely an A-star. iKilic et al.l ( 120091 . 12010a ). iKawka et al.l ( 120101 ). and ISteinfadt et al. 



( 120101 ) describe radial velocity observations of two more ELM WDs identified in the litera- 
ture. Thus, there are 23 (candidate) ELM WD systems with comprehensive radial velocity 
measurements. All but three of these systems (J1448+1342, J2252-0056, and J2345-0102) 
are in binary sytems. 

Limiting this sample to 19 WDs with M < 0.25M Q , 18 are in short period binary 
systems. Thus, close binary evolution is required for the formation of ELM WDs. The only 
system that does not display significant radial velocity variations, J1448+1342, is at the 
upper mass limit of this selection; all 18 WDs with M < O.25M are indeed binaries. The 
average velocity semi-amplitude of these 18 binaries is 236 km s _1 , whereas the upper limit 
for the velocity semi-amplitude of J1448+1342 is 35 km s _1 . An average system viewed 
at i < 8.5° would be consistent with the observations of J1448+1342. For a randomly 
distributed sample of orbital inclinations, there is a 1.1% chance that i < 8.5°. Thus, there 
is a 21% likelihood of finding one of the 19 sys tems with i < 8.5° . It is possible that 



J1448+1342 is a pole-on binary system (see also IBrown et al.l l2010al ). Hence, the binary 
fraction of M < 0.25M© WDs is at least 95% and it may be as high as 100%. 



Th is result is in contr a st to t he bin ary fraction o f slight ly more massive WDs. iMaxted et al 



feOOOl ). iNapiwotzki et al.l ( 120041 ) . and IBrown et al 
of M w O.4M WDs are single (see also IKilic et al. 



(2010) 



2010c 



demonstrate that about 30-50% 
Si ngle low-mass WD s may form 



through enhanced mass loss from a metal- rich progenitor star ( iKilic et al.ll2007d ). but we do 
not expect even the most metal-rich stars to loose enough mass to create ~ O.2M WDs. 
Here we confirm this expectation that the lowest mass WDs all form in close binary systems. 



4.2. A Dozen Merger Systems 



Double degenerate (DD) mer ger systems are one of the two main accep t ed formation 



channels for Ty pe la supernovae (llben fc Tutukovl 1 19841 ; iBranch et al.l Il995l ; iKotakl 12008 



Di Stefanoll2010l . and references therein). However, only a handful of DD systems are known 
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to have combined mass above the Chandrasekhar limit and short enough orbital per i ods to 
merge through gravitational wave radiation within a Hubble time. iNapiwotzki et al.l (120011 ) 
started the heroic SPY survey for short period DD systems. They used the Very Large 
Telescope to observe 1014 stars brighter than 16.5 mag, in c ludin g 75% of the known WDs 
accessible from the southern hemisphere. INapiwotzki et al.l (120041 ) report that the SPY sur- 
vey discovered a bout 100 new DE > syste m s, one of which is a potential SNe la progenitor 
(WD 2020-425. INapiwotzki et al.ll2007h . beier et al.l fj201oh describe the discovery of an- 
other massive DD merger system and also a summary of the current discoveries from the 
SPY survey (their Figure 4). 

Figure 10 shows the total masses and periods for the short period DD systems (including 
sdB+WD) discovered by the SPY survey (kindly made available to us by S. Geier) and 
our ELM WD sample. To explain the observed SNe la rate, only one out of 1000 WDs is 
expected to be a SNe la progenitor (INelemans et al.ll200ll ). Hence, the relatively few number 
of potential SNe la progenitors in the SPY survey is not surprising. What is striking is that 
the number of merger systems found in the SPY survey is significantly smaller than our ELM 
WD survey. We discovered 12 binary WD merger systems (r morge < r hubble ) in a sample of 
23 stars, whereas the SPY survey discovered a handful of such systems in a sample of 1000 
stars. The reason for this is the binary common envelope evolution required to form ELM 
WDs. 



Figure 11 shows the total masses and periods for previously identi fied WD+WD binary 
systems in the literature (INelemans et al.ll2005b INapiwotzki et al.ll2007l ). There were only six 
double WD merger systems known prior to our work. With the additional 12 merger systems 
presented in Table 4, we now have tripled the number of WD merger systems known. 



4.3. The Period Distribution of Binary WDs 



Nelemans et al.l ( 1200 ll ) studied the period and mass distribution of double WD systems 



using population synthesis models. The details of these models including selection effects 
due to orbital evolution of the systems (disappearence of the shortest period systems due 
to mergers), target select ion biases in the curren t surveys, and different cooling models are 
discussed in their paper. INelemans et al.l ( 120011 ) find a correlation between the masses and 
periods; more massive WD primaries have longer periods. This is not surprising; shorter 
period systems would start interacting earlier in their evolution compared to longer period 
systems and experience enhanced mass-loss during the red giant phase, hence end up as 
lower mass WDs. They also find that the mass ratio distribution peaks at q — 1, favoring 
equal mass binary components. This distribution is consistent with the observed mass ratio 
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distribution of binary WDs known prior to our work. 

Figure 12 shows the masses and orbital periods for our sample of ELM WDs and the 
previously kn own binary WD s ystems in th e literature compared to the population synthe- 
sis models of iNelemans et al.l ( 120011 . 120041 . kindly made available to us by G. Nelemans). 
Depending on which cooling model is used, the model predictions change somewhat (see 
Nelemans et al.l 12004 for a discussion of the cooling models used here). However, the main 
trend in the mass versus period distribution is evident; lower mass primaries are in shorter 
period systems. This result is consistent with the observed distribution of masses and pe- 
riods for previously known binary WD systems and the ELM WD systems. All but one of 
the M < O.3M systems have periods shorter than a dajfl, whereas the majority of the more 
massive WDs have periods longer than a day. 

The overall number distribution of the population synthesis models and the observations 
cannot be directly compared without taking into account the obvious biases in the target 
selection. The sample of previously known binary WDs are selected from a magnitude 
limited sample down to 16-17 mag, whereas our sample mainly comes from fainter targets 
selected from the SDSS, which has its own spectroscopic target selection bias. Nevertheless, 
an overabundance of 0.17 M W Ds is ev i dent in Figure 12. A comparison between the 



popul ation synthesis models using iHansenl (119991 top panel) and the modified iDriebe et al. 



( 119981 . bottom panel) cooling models shows that the expected number of the lowest mass 
WDs changes depending on the cooling model u sed. The latter mod els predict a significantly 
larger number of 0.17 M WDs. Based on the IPanei et all (120071 ) models, a 0.25 M WD 
cools 4x faster than a 0.17 M WD down to an absolute magnitude of M g = 10 mag. Hence, 
the observed mass distribution of ELM WDs in a magnitude-limited survey is skewed toward 
lower mass objects (IBrown et al.ll2010bl ). Population synthesis mo dels and a large, unbi ased, 
magnitude-limited survey of low-mass and ELM WDs like that of IBrown et al.l ( j2010al ) will 
be useful to constrain the ELM WD cooling models. 

Another important prediction of the population synthesis models is that there should 
not be many systems with periods less than an homo We have three ELM WD systems 
with periods close to 1 hr. The shortest period system, J1053+5200 has a period of 61 
min ( iKilic et al.ll2010bl ; iMullally et al.ll2009l ). The shorte st period hydrogen-rich cataclysmic 
variable (CV) V485 Cen has an orbital period of 59 min (lAugusteijn et al.lll996l ). Both short 
period binary WDs and CVs are expected to contribute to the AM CVn population, which 



4 The single velocity non-variable system, J1448+1342, is probably a pole-on or a relatively long-period 
binary system. 

5 These systems either merge or evolve into AM CVn systems relatively quickly. 
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contains interacting WD binary systems with orbital periods ranging from a few minutes to 
about an hour. The two longest per iod AM CVn systems have periods of 56 min and 66 
min, respectively (IRoelofs et al.ll2007l ). It seems that hydrogen- deficient AM CVn interacting 
binary systems have periods less than about an hour and the hydrogen-rich detached WD 
systems have orbital periods about an hour or longer. 



4.4. The Minimum Mass for Low-mass WDs 



Before the SDSS, ELM WDs were traditionally found as companions to MSPs. The 
masses for the MSP companions are usually measured from the pulsar ti ming measurements 
since the WDs are usual ly too faint for follow-up optical spectroscopy. Ivan Kerkwijk et al. 



(119961 ) and lBassa et al.l (120061) obtained Keck and VLT spectroscopy of companions to two 
MSP systems. Ivan Kerkwiik et~aD f)l996h derive T pff = 8550 ± 25 K and log g = 6.75 ± 0.07 
for the companion to PSR J1012+5307, whereas iBassa et al.l (120061 ) obtain T e g- = 10,090 ± 



150 K and log g = 6.44 ± 0.20 for the companion to PSR J1911-5958A. These two WDs are 
included in Figure 8. 



Kilic et al.l ( l2007aU bl) discov ered the lowes t surfac e gravity WD known at the time, SDSS 



J0917+4638, with log g = 5.55. iBrown et al.l ( 12010al ) discovered two WDs with even lower 



( 


2001, 


2002 


) and 


Panei et al. 



(120071 ) models, all of these WDs, including the two MSP companions, have M > O.17M . 
There are about a dozen WDs in Figure 8 with M « O.17M , but none below that mass 
limit. However, this result is most likely due to observational biases present in the SDSS and 
the Hypervelocity Star Survey. The former is, of course, an incomplete spectroscopic survey 
and the latter h as a color selection t hat excludes T U < 10000 K objec ts with log g = 5 — 6 
(see Figure 1 in IBrown et al.ll2010al ). Based on the iPanei et al.l (120071 ) models, a 0.16 M 



WD takes 15 Gyr to cool down from 9200 K to about 7000 K with an increase in surface 
gravity from log g = 5.7 to 6.4 and a decrease in luminosity from M g = 8.1 to 11.1 mag. 
Such WDs are not included in the Hypervelocity Star Survey color selection. 

Based on radio timing measurements, 45 out of 99 MSPs with periods less than 20 ms 
have companions less massive than O.17M (B. Kiziltan 2010, private communication). The 
spectral types of the companions are unknown. However, a large fraction of these objects are 
expected to be low- mass WDs; nature produces M < 0.1 7M WDs in MSP systems. Lower 
mass ELM WDs should be produced in binary WD systems as well. Population synthesis 
calculations for double WDs do not predict many WDs below 0.17 M , but such WDs may 
eventually be found in double WD systems in the SDSS or other spectroscopic surveys. 
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4.5. Merger Products: (underluminous) SNe la, AM CVn, and RCrB stars 



The end result of the future evolution of the 12 ELM WD merger systems in our sample 
depends on the mass ratios of the two components. The systems with extreme mass ratios 
(q < 0.2) should have stable mass transfer at initial Roche lobe filling. Stable mass transfer 
leads to the formation of AM CVn systems that transfer angular momentum back to the 
orbit, which increases the orbital period. If the higher mass WDs in such systems have masses 
close to the Chandrasekhar limit, accretion may result in a Type la SNe explosion. However, 
the probability of such an event is small (< 7%) for our systems. This probability would 
be higher if sub- Ch andrasekhar ma s s WD explosions al so contribute to the observed Type 
la SNe population (ISim et al.ll2010l ; iKromer et al.l |2010L and references therein). However, 
our ELM survey is biased against SNe la progenitors due to the selection of very low-mass 
primaries. 

Figure 13 shows the mass ratios and periods for our ELM WD sample plus the previously 
known binary WD systems. Single-lined systems are shown as filled triangles, whereas the 
double-lined systems and the eclipsing WD system NLTT 11748 are shown as filled circles. 
Open circles mark the objects with merger times shorter than a Hubble time. Out of the 18 
merger systems, six have mass ratios smaller than 0.25 (assuming an average inclination of 
60°), and are likely to form AM CVn systems. These six systems also have the highest total 
mass (M1+M2) in our sample. J1233+1602 is the most extreme case with Mi+M 2 = 1.37M Q 
for i = 60°. It is striking that the number of extreme mass ratio systems in our sample (six) 
i s the same as the number of AM CVn systems found in the SPS S Data Relase 5 are a 
(IRoelofs et al.ll2007l ). These systems are potential SNe .la progenitors (IBildsten et al.ll2007l ). 



Based on a magnitude-limited sample of six ELM WD mergers, iBrown et al.l (l2010bl ) 
find that half of their sample has mass ratios of r^O.15 and they estimate a contribution 
factor of >3% to the AM CVn formation rate from the observed population of ELM WDs. 
The fraction of extreme mass ratio systems i n our survey i s also 50% (six out of 12 ELM 
WDs), supporting the conclusions reached by IBrown et al.l (j2010bl ). 



The o t her six ELM WD systems in our sample have mass ratios larger than 0.4. 
Kilic et al.l (j2010al ) discuss various evolutionary scenarios for these objects. These systems 
will merge within a Hubble time and form a variety of exotic objects including underlumi- 
nous SNe, extreme helium stars (RCrB), or single helium-enriched subdwarf O stars. The 
merger rate of ELM WDs is roughly the same as the observed rate of underluminous SNe 
farown et alihoiObh . 



SDSS J1234— 0228 is a unique ELM WD merger system. There is a 85% chance that the 
combined mass of the two stars in this system is lower than 0.45 M & . Therefore, within the 
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next few billion years, this system will merge and likely form a a single low-mass WD. There 
are several scenarios to explain the presence of si ngle low-mass (~ O.4M ) WDs, including 
enhanced mass-loss from metal-rich progenitors f lKilic et al.ll2007d ) and mass loss through 
common envelope evolutio n with a low-mass brown dwarf or planetary mass companion 
( iNelemans fc Taurid Il998l). Single low- mass WDs can also form through mergers of two 
lower-mass WDs (jlben et al.l Il997t ISaio fc Jefferyl 120001 ). though no such progenitors have 
been found until now. Here we have uncovered a potential progenitor system for a future 
single low-mass WD, J1234-0228. 



About 10% of DA WDs are low-m ass ( ILiebert et al.l 120051 ) and half of these low- mass 
WDs are single (Napiwotzki et al.ll2004[ ). corresponding to a space density of 2.5 xlO 5 kpc~ 3 
rtSion et al.ll2009h. The SPS S DR4 WD catalog covers 4783 square degrees with g = 15 - 20 
mag (jEisenstein et al.l 120061 ). We estimate the loca l space density using the modified l/V max 
method and the disk model of iRoelofs et al.l (120071 ). Based on one system, J1234— 0228, the 
space density of merging WD systems that can form single low- mass WDs is 1.1 kpc~ 3 . This 
number is likely uncertain by a factor of 10 due to incompleteness of the SDSS spectroscopic 
target selection. Nevertheless, it is several orders of magnitude smaller than the estimated 
space density of single low-mass WDs. Hence, this formation channel does not significantly 
contribute to the single low-mass WD population in the solar neighborhood. 



5. CONCLUSIONS 

Our MMT ELM survey has now discovered 20 short period binary WD systems, in- 
cluding 12 merger systems. Only six binary WD merger systems were known prior to our 
work. The ELM survey has now tripled the number of such systems. These merger systems 
provide potential formation channels for a variety of interesting objects including underlu- 
minous SNe, AM CVn, extreme helium stars, single subdwarf stars, and now single low-mass 
WDs. 

All but one of the 19 stars with M < O.25M are in binaries, strongly arguing for a 
binary formation channel for ELM WDs. The combined sample of more than 40 short period 
WD systems found in our ELM survey and the literature demonstrate that the shortest 
period detached WD systems have ~1 hr orbital periods and they also have the lowest mass 
primaries. In addition, the lowest mass WDs in double WD systems have M = O.17M . 
Lower-mass objects may be produced in such systems but are currently undetected in the 
SDSS and other spectroscopic surveys. 

The sample of ELM WD merger systems will be significantly enlarged by our continuing 
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ELM survey based on SDSS photometry. Such a magnitude-limited survey will not suffer 
from the biases present in our current sample, and will provide invaluable information on 
the merger rates of ELM WD systems and their contribution to underluminous SNe and AM 
CVn systems. Furthermore, a large sample of ELM WDs and population synthesis models 
will improve our understanding of the cooling models for ELM WDs. 

We thank G. Nelemans for providing his binary WD population synthesis model results, 
S. Geier for the latest results from the SPY survey, and D. Koester for WD model spectra. 
Support for this work was provided by NASA through the Spitzer Space Telescope Fellowship 
Program, under an award from Caltech. 

Facilities: MMT (Blue Channel Spectrograph), CXO (ACIS-S) 
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Table I. X-ray Observations of ELM WDs 



Name 


ObsID 


Expos. 


Dist. 


N H 


Countrate 


Lx 


(SDSS) 




(ks) 


(kpc) 


(cm- 2 ) 


(cts s" 1 ) 


(ergs s _1 ) 


J105353. 89+520031.0 


9963 


20.50 


1.08 


1.0 x 10 20 


< 2.2 x 10~ 4 


< 9.5 x 10 28 


J123410.36-022802.8 


9964 


5.50 


0.78 


2.2 x 10 20 


< 8.4 x 10~ 4 


< 2.0 x 10 29 


J092345. 59+302805.0 


9965 


1.03 


0.27 


2.0 x 10 20 


< 4.5 x 10~ 3 


< 1.3 x 10 29 



Note. — 99% confidence X-ray count rate upper limits for three ELM WDs from Chandra 
X-ray observations. Countrate limit is in 0.3-6 keV band and Lx limit is in 0.5-6 keV band. 



Table 2. Orbital Parameters 



Object 


P 


K 


^systemic 


Spec. Conjunction 


Mass Function 




(days) 


(km s- 1 ) 


(km s- 1 ) 


(HJD - 2454000) 


(M ) 


J0022+0031 


0.4914 ± 0.0254 


80.8 ± 1.3 


-20.3 ± 0.8 


732.71818 


0.02681 ± 0.00190 


J0022-1014 


0.0799 ± 0.0030 


145.6 ± 5.6 


-38.5 ± 3.7 


732.76530 


0.02553 ± 0.00310 


J1234-0228 


0.0914 ± 0.0040 


94.0 ± 2.3 


+50.3 ± 1.8 


917.75366 


0.00787 ± 0.00067 


J1625+3632 


0.2324 ± 0.0396 


58.4 ± 2.7 


-95.0 ± 2.1 


922.74700 


0.00480 ± 0.00105 


J2252-0056 






-23.2 ± 0.8 






J2345-0102 






-161.2 ± 0.4 







Note. — J0022+0031 has a period alias at 7.9 hr. 



Table 3. Physical Parameters 



Object 


T e ff 


log g 


Mass 


M 2 


M 2 (i = 60°) 


NS 


SN la 


Tmergc 




(K) 




(M Q ) 


(M ) 


(M G ) 


Prob. 


Prob. 


(Gyr) 


J0022+0031 


17890 ± 110 


7.38 ± 0.02 


0.38 


> 0.21 


0.26 


2.5% 


1.9% 


< 75.7 


J0022-1014 


18980 ± 380 


7.15 ± 0.04 


0.33 


> 0.19 


0.23 


2.3% 


1.4% 


< 0.73 


J1234-0228 


18000 ± 170 


6.64 ± 0.03 


0.23 


> 0.09 


0.11 


0.9% 


0.3% 


< 2.69 


J1625+3632 


23570 ± 440 


6.12 ± 0.03 


0.20 


> 0.07 


0.08 


0.6% 


0.2% 


< 45.2 


J2252-0056 


19450 ± 270 


7.00 ± 0.02 


0.31 












J2345-0102 


33130 ± 450 


7.20 ± 0.04 


0.42 













Table 4. Merger Systems in the ELM Survey 



Object 


T e ff 


logs 


P 


K 


Mass 


M 2 


M 2 (i = 60°) 


^mcrgc 




(K) 




(days) 


(km s -1 ) 


(M e ) 


(M G ) 


(Me) 


(Gyr) 


J0022-1014 


18980 


7.15 


0.07989 


145.6 


0.33 


> 0.19 


0.23 


< 0.73 


J0755+4906 


13160 


5.84 


0.06302 


438.0 


0.17 


> 0.81 


1.12 


< 0.22 


J0818+3536 


10620 


5.69 


0.18315 


170.0 


0.17 


> 0.26 


0.33 


< 8.89 


J0822+2753 


8880 


6.44 


0.24400 


271.1 


0.17 


> 0.76 


1.05 


< 8.42 


J0849+0445 


10290 


6.23 


0.07870 


366.9 


0.17 


> 0.64 


0.88 


< 0.47 


J0923+3028 


18350 


6.63 


0.04495 


296.0 


0.23 


> 0.34 


0.44 


< 0.13 


J1053+5200 


15180 


6.55 


0.04256 


264.0 


0.20 


> 0.26 


0.33 


< 0.16 


J1233+1602 


10920 


5.12 


0.15090 


336.0 


0.17 


> 0.86 


1.20 


< 2.14 


J1234-0228 


18000 


6.64 


0.09143 


94.0 


0.23 


> 0.09 


0.11 


< 2.69 


J1436+5010 


16550 


6.69 


0.04580 


347.4 


0.24 


> 0.46 


0.60 


< 0.10 


J2119-0018 


10360 


5.36 


0.08677 


383.0 


0.17 


> 0.75 


1.04 


< 0.54 


NLTT 11748 


8690 


6.54 


0.23503 


273.4 


0.18 


0.76 




7.20 



Note. 
d2010ah. 



and 



Kilic et 



This table is based on the data p resented in Kilic et al 
d2010ch . and this paper. See lMullallv et al.l |2009T T 



al 



Kawka et al 



( 2010bl), Brown et al 



Steinfadt et al.l (|2010h . 



(20101) for additional observations of J1053+5200, J1436+5010, and NLTT 



11748. 
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Fig. 1. — Radial velocities and the best-fit orbit for J0022+0031. The bottom right panel 
shows all of the data points phased with the best-fit period. 
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Fig. 2.— Same as Figure 1, but for J0022-1014. 
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Fig. 3.— Same as Figure 1, but for J1234-0228. 
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Fig. 4.— Same as Figure 1, but for J1625+3632. 
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Fig. 5.— Same as Figure 1, but for J2252-0056. 
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Fig. 6.— Same as Figure 1, but for J2345-0102. 
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Fig. 7. — Model fits (red solid lines) to the composite spectra of our targets (jagged lines). 
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Fig. 8. — The best fit solutions for the surface gravit y and tempe r ature s of our targets (red 
filled circles) , overl aid on tracks of constant mass from lPanei et al.l (120071 ) models updated by 
Kilic et al.l ( I2010a[ ). The dashed and dotted lines show solar metallicity and halo metallicity 
(Z=0.001) models of ISerenelli et al.l toOli I2OO2I ) for 0.17 M WDs, respective ly. Spectro- 



scopi cally confirmed WDs from the ELM survey and the sdB star HD 188112 ( iHeber et al. 



20031 ) are shown as black circles. Companions to milli-second pulsars PSR J1012+5307 and 
J1911-5 958A are shown a s triangles. The four objects with X-ray nondetections (Table 1 
and lAgiieros et al.ll2009bl ) are marked by open circles. 
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Fig. 9. — Best-fit model (red solid line) to the observed spectrum of SDSS J2049+0005. The 
best-fit model has log g = 5 and does not match the higher order Balmer lines. J2049+0005 
is likely an A star. 
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Fig. 10. — Total mass plotted again st logarithmic period of double degenerate systems from 
the SPY survey (IGeier et al.l 120101 ) and our ELM WD sample (red filled circles). ELM 
WDs, central stars of planetary nebulae (CSPN), double- lined /single- l ined W D systems, 



and sdB/sdO systems are shown with different symbols. IGeier et al.l ( 120101 ) assumes an 
inclination angle of 52° for the single-lined systems, whereas we use an inclination angle of 
60° for our systems. 
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Fig. 11. — Same as Figure 10, but including only the double WD systems found in the 
literature (INelemans et al.ll2005l ; iNapiwotzki et al.l 120071 ). 
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Fig. 12. — Model population of double WDs as a function of orbital period and mass of 
the brighter W D of the pair (kindly made available by G. Nelemans). The top panel uses 
Hansen! (119991 ) cooling models whereas the bottom pa nel uses modified iDriebe et al.l (119981 ) 
cooling models. The observed binary WD population (INelemans et al.ll2005l ) and our ELM 
WD sample (red filled circles) are also shown. 
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Fig. 13. — Mass ratios and periods for binary WD systems found in our ELM survey and 
the literature. The mass ratio limit for stable mass transfer is also indicated. 
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APPENDIX A. Radial Velocity Data 
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Table 5. Radial Velocity Measurements 



Object 


HJD 


Vhelio 




-2450000 


/i _i \ 

(km s i ) 


T0022+0031 


732 722684 


—72 2 + SO ^ 




732 85851 Q 






732 Q3QQ55 


91 9 4- 36 




732 Q574RR 


— 1 Q 3 + 50 4 




733 72825Q 


_ak 4 x in t; 




733 85023Q 


1 C O _|_ 1 1 fi 




733 Q51 3Q7 


-1(19(1 + 47Q 




1 0Q4 846527 


— 12.4 ± 12.4 




1 100 810Q77 

J-J-lyly.O-LVJiyf 1 


en 7 _|_ ICC 




1 100 827123 


•{0 7+ 1Q fi 




1 100 841278 






1 1 00 858385 


sn fi + 14 s 




1 1 00 872447 

_1_ J_ Vy Vy . O 1 Z/^t4: 1 


74 fi _i_ i q e 

1 ^.U -I— -L t/ .O 




1 1 00 8Q0352 


q<> c _i_ ico 




1 1 00 Q04427 

J_ J_ Vy Vy . i/Ut:t:Zi 1 






1 100 Q1Q5fifi 


51 q 4. oq 7 




1 1 01 8Q545Q 


i q q 4. 99 q 




1 1 01 Q0R050 


Qfj q 4. 91 « 




1 1 01 Q1 778fi 

J- 1U J- . Cf J- i 1 ou 


J. 1 . U _l_ Z( _I_ . Vy 




1 1 02 Q04Q28 


qo 9 4. 1 q q 




1 1 02 Q1 5530 


wi 4- 9fi q 




1 384 044307 


1 nn fi _i_ 90 7 






— < 0.1 ± 41). 




1387.925074 


-14.8 ± 72.2 




1395.943892 


-96.8 ± 19.8 




1396.936261 


-87.0 ± 17.1 




1396.976913 


-11.9 ± 49.8 


J0022-1014 


732.748459 


-156.9 ± 27.6 




732.759199 


-59.1 ± 41.4 




732.832753 


-232.5 ± 32.4 




732.843424 


-77.3 ± 32.2 




732.913470 


-199.3 ± 24.1 




732.924142 


13.5 ± 20.7 




733.781989 


-198.3 ± 44.7 




733.872718 


-125.5 ± 47.9 
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Table 5 — Continued 



Object 


HJD 


uhelio 




-2450000 


(km s- 1 ) 




1094.862244 


25.5 ± 14.1 




1102.756314 


89.0 ± 31.1 




1102.770423 


67.9 ± 31.6 




1102.787923 


-148.8 ± 20.6 




1102.801997 


-214.7 ± 22.9 




1102.818895 


33.5 ± 18.9 




1102.832969 


59.5 ± 15.6 




1102.848212 


46.9 ± 18.2 




1102.862274 


-121.9 ± 19.0 




1102.876371 


-148.9 ± 33.8 


J1224-0228 


917.805724 


-1.6 ± 6.2 




922.728604 


98.1 ± 6.4 




922.815478 


111.6 ± 5.8 




922.904031 


132.2 ± 14.9 




923.754131 


-16.3 ± 5.2 




923.867625 


-22.5 ± 4.6 




923.931236 


6.6 ± 7.9 




924.850731 


52.6 ± 13.4 




924.857965 


-61.9 ± 17.8 




1275.829077 


136.3 ± 10.1 




1275.834841 


73.9 ± 11.9 




1275.841936 


17.8 ± 10.5 




1275.847700 


-20.6 ± 7.8 




1275.853441 


-15.4 ± 9.3 




1275.860270 


-13.8 ± 11.7 




1275.866022 


-21.6 ± 11.9 




1275.871763 


-27.3 ± 8.4 




1275.919252 


146.8 ± 11.0 




1275.925560 


103.6 ± 12.5 


J1625+3632 


922.872132 


-90.0 ± 12.8 




922.957182 


-146.1 ± 9.7 




923.002704 


-39.0 ± 9.6 




923.837051 


-160.2 ± 10.4 




923.912980 


-90.8 ± 8.4 




923.956569 


-54.9 ± 7.6 




924.004812 


-73.7 ± 16.0 
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Table 5 — Continued 



Object HJD 


uhelio 


-2450000 


(km s- 1 ) 


925.000912 


-134.7 ± 11.0 


925.012835 


-140.7 ± 21.9 


1276.936718 


-12.2 ± 14.2 


1276.943038 


-25.0 ± 17.5 


1276.951372 


-11.3 ± 17.6 


1276.957668 


-57.6 ± 17.7 


1276.964011 


-60.2 ± 13.3 


1276.972808 


-7.3 ± 24.7 


1276.979128 


-34.2 ± 27.0 


1276.985447 


-50.1 ± 26.9 


1384.783927 


-51.8 ± 20.6 


1385.738969 


-59.6 ± 11.9 


1396.823036 


-23.7 ± 42.7 


J2252-0056 732.661721 


-33.0 ± 8.4 


732.706500 


-52.3 ± 13.5 


732.786660 


-40.0 ± 11.8 


732.872943 


-17.0 ± 16.2 


733.653315 


-24.9 ± 16.9 


733.752698 


4.8 ± 13.1 


733.761692 


-39.8 ± 13.9 


733.829687 


-18.2 ± 18.7 


733.892232 


-11.5 ± 11.7 


1094.801506 


-22.9 ± 7.3 


1100.666435 


-18.0 ± 16.4 


1100.681967 


4.7 ± 9.2 


1100.696261 


-14.5 ± 8.9 


1100.712232 


34.4 ± 8.3 


1100.726399 


-12.9 ± 9.5 


1100.742289 


1.4 ± 8.7 


1100.756386 


26.4 ± 10.8 


1100.770645 


-29.4 ± 15.4 


1100.781270 


-29.9 ± 11.5 


1101.650617 


-35.3 ± 15.1 


1101.659471 


20.6 ± 13.9 


1101.672399 


-118.1 ± 16.3 


1101.682989 


27.1 ± 16.7 
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Table 5 — Continued 



Object HJD 


Vhelio 


-2450000 


(km s _1 ) 


1101.698706 


-35.8 ± 13.6 


1101.709307 


-2.8 ± 16.8 


1101.719897 


-33.6 ± 17.3 


1102.645214 


-86.1 ± 19.9 


1102.655804 


28.9 ± 20.7 


1102.667655 


-65.7 ± 23.2 


1102.678245 


-63.8 ± 15.3 


J2345-0102 732.689649 


-219.3 ± 27.1 


732.799579 


-130.2 ± 38.3 


732.810239 


-94.4 ± 35.0 


732.889312 


-114.6 ± 30.3 


733.697952 


-129.3 ± 31.9 


733.711030 


-159.5 ± 53.8 


733.809108 


-142.7 ± 24.4 


733.928533 


-153.0 ± 47.3 


1094.830335 


-175.5 ± 28.8 


1101.735791 


-139.1 ± 59.8 


1101.751161 


-242.8 ± 64.4 


1101.765224 


-243.8 ± 41.6 


1101.786080 


-109.0 ± 31.3 


1101.800154 


-97.8 ± 23.5 


1101.816149 


-199.4 ± 40.1 


1101.830211 


-138.7 ± 25.3 


1101.846137 


-109.9 ± 26.4 


1101.860199 


-138.8 ± 22.3 


1102.701680 


-136.8 ± 29.2 


1102.719307 


-95.9 ± 65.7 


1102.736402 


-135.1 ± 46.7 



